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Abstract:  This  paper  presents  a  general  existence  and  uniqueness  theory  for  differential-algebraic  equations  ex¬ 
tending  the  well-known  ODE  theory.  Both  local  and  global  aspects  are  considered,  and  the  definition  of  the  in¬ 
dex  for  nonlinear  problems  is  elucidated.  For  the  case  of  linear  problems  with  constant  coefficients  the  results 
are  shown  to  provide  an  alternate  treatment  equivalent  to  the  standard  approach  in  terms  of  matrix  pencils. 
Also,  it  is  proved  that  general  differential-algebraic  equations  carry  a  geometric  content,  in  that  they  are  locally 
equivalent  to  ODEs  on  a  *  constraint"  manifold.  A  simple  example  from  particle  dynamics  is  given  to  illustrate 
our  approach.  0 


1.  Introduction 


Differential -algebraic  equations  (DAEs)  are  frequently  identified  as  implicit  equations 


F(t,x,x')  =  0  (1.1) 

for  which  the  derivative  x  cannot  be  expressed  explicitly  as  a  function  of  t  and  x  (see  e.g.  [1]).  In  particular,  if 
x  e  R“  and  F  maps  into  R",  this  includes  the  case  when  the  partial  derivative  DpF(t , x, p)  of  F  with  respect 
to  its  third  variable  p  is  not  surjective.  More  specifically,  in  the  setting  of  DAEs  it  is  natural  to  require  the 
stronger  hypothesis  that  DpF(t,  x,  p)  has  constant  rank  on  the  domain  under  consideration.  Indeed,  the  proto¬ 
type  for  such  equations  is  given  by 


F(t,x,p)  = 


F2(t,x,p) 


(12) 


where  Fj  and  F%  map  into  R'*-''  and  R'',  respectively,  and  DpF2(t,  x,  p)  has  full  rank,  so  that,  indeed, 
DpF(t,  x,  p)  has  constant  rank  r  <  n. 

Many  DAE-problems  of  practical  interest  do  not  exhibit  such  a  convenient  splitting  between  algebraic  and 
differential  parts  as  in  (1.2).  Moreover,  even  if  the  equations  can  be  whiten  in  the  separated  form  (1.2),  the  rank 
of  DpF2(t ,  x,  p)  may  turn  out  to  be  less  than  r  so  that  F2(t,  x,  x')  =  0  is  an  equation  containing  an  implicit 


1)  Thu  work  waa  supported  in  pan  by  O NR -gram  N -00014-90-1-1025  and  NSF-granl  CCR-8907654. 
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algebraic  part.  While  these  comments  suggest  the  need  for  a  thorough  investigation  of  DAEs  in  the  broad  setting 
of  (1.1),  existence  and  uniqueness  theories  for  these  equations  have  not  been  developed  in  such  generality  and 
available  results  place  more  restrictive  conditions  on  the  form  of  F  (see  e.g.  [4],  [7], [8]). 

In  this  paper  we  present  existence  and  uniqueness  results  for  general  problems  of  the  form  (1.1)  nmw 
some  "generic"  conditions  which  do  not  assume  that  F  has  a  special  form.  Briefly,  our  technique  consists  in 
deriving  a  necessary  condition  for  the  existence  of  local  solutions  by  differentiating  the  equation  once  followed 
by  an  application  of  orthogonal  projections  onto  the  range  of  DpF  to  produce  a  system  that  again  contains  only 
x  and  x  .  For  this  decomposed  system  a  sufficient  condition  is  then  obtained  that  guarantees  the  equations  to 
define  a  local  vector  field  to  which  the  standard  existence  theory  applies.  The  effectiveness  of  this  conceptually 
-  but  not  always  technically  --  simple  approach  for  providing  an  answer  to  this  notoriously  complex  problem 
may  be  called  surprising. 

The  sufficient  condition  essentially  requires  that  the  index  (see  e.g.[l])  of  the  DAE  is  one.  But  it  also 
turns  out  that  the  theory  can  be  applied  recursively  to  equations  of  higher  index  provided  that  the  resulting  equa¬ 
tions  again  satisfy  the  constant  rank  condition  mentioned  above.  This  is  a  requirement  of  a  global  nature  not 
covered  by  the  standard  index  theories.  In  fact,  it  also  suggests  that  the  constant  rank  condition  is  inherent  to  the 
definition  of  the  index. 

In  Section  2  below  we  present  the  local  existence  and  uniqueness  theory  sketched  above.  Then  in  Section 
3  these  local  solutions  arc  extended  under  conditions  which  correspond  to  those  of  the  standard  ODE  theory.  In 
Section  4  it  is  proved  that  when  DF  has  full  rank  then  the  set  of  admissible  initial  points  forms  an  r- 
dimensional  submanifold  of  R'xR"  and  the  DAE  is  locally  equivalent  to  a  differential  equation  on  an  r- 
dimensional  submanifold  of  R".  This  also  shows  that  the  geometric  approach  developed  in  [7]  and  [8]  is  con¬ 
ceptually  valid  in  general.  Then  in  Section  S  we  apply  our  results  to  linear  equations  with  constant  coefficients 
and  prove  that  the  recursive  application  of  the  technique  leads,  exactly  to  the  standard  index  for  such  linear 
DAEs.  Finally,  Section  6  concerns  the  generalization  of  this  recursive  application  of  the  results  to  the  general 
nonlinear  case  which,  as  mentioned  before,  requires  the  additional  global  assumption  that  the  constant-rank  con¬ 
dition  remains  valid.  This  is  illustrated  on  the  classical  example  of  the  nonlinear  pendulum. 
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2.  Local  Existence  and  Uniqueness 

For  case  of  notation,  wc  shall  consider  (1.1)  first  in  the  autonomous  form 

F(x,x')  =  0  (2.1) 

where  it  is  assumed  that 

f:  R*  xR"  - »  R"  is  of  class  C2  on  the  open  set£:  c  R"  x  R*  (2.2a) 

and 

rank  DpF(x,  p)  =  r  <  n ,  for  all  (x ,  p )  e  £.  (2.2b) 

Our  results  will  show  that  the  differentiability  assumption  (2.2a)  --  instead  of  the  expected  and  apparently 
more  natural,  minimal  C ‘-regularity  of  F  —  turns  out  to  be  important  for  the  theory.  We  note  also  that  the  char¬ 
acter  of  the  problem  changes  significandy  when,  instead  of  (2.2b),  DpF(x,  p)  is  only  required  to  be  singular  on 
some  lower  dimensional  sub-manifold  of  R  *  x  R  *  (sec  [5]).  Finally,  we  observe  that,  while  the  rank  condition 
(2.2b)  may  suggest  a  transformation  of  the  problem  to  some  canonical  form  by  means  of  a  version  of  the  "rank 
theorem"  (see  [3]),  it  must  be  noted  that  such  a  reduction  mixes  the  variables  x  and  p.  Since  in  (2.1)  the  deriva¬ 
tive  x’  occupies  the  position  of  the  p-variable,  such  a  reduction  is  not  readily  usable  to  transform  (2.1)  to  an 
explicit  ordinary  differential  equation. 

A  C2-solution  of  (2  1)  shall  be  any  function 

x:  J  R"  ,  (x(t),  x'(t))  e  £,  for  t  e  J,  (23) 

which  is  of  class  C2  on  some  open  interval  J  cR"  and  satisfies  F(x(t),  x'  (f))  =  0  for  all  t  e  /.  For  any  Ca- 

solution  (2.3)  of  (2.1)  wc  obtain  by  differentiation 

>  ■  »  >  »  ■  •  » 

DxF(x(t),  x' (t))x'(<)  +  DpF{x(t),  x7 (t))x"(t)  =  0,  t  e  J,  (2.4) 

which  provides  the  following  necessary  condition: 

Lemma  2.1:  If  (2.2a)  holds  for  (2.1)  then  for  a  given  point  (x,  p)  e  £  the  two  conditions 

£(x,p)  =  0  (2.5a) 


DtF(x,p)p  e  rge  DpF(x,p) 


(2.5b) 
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are  necessary  for  the  existence  of  a  C:-soluuon  (2.3)  of  (2.1)  that  passes  through  (x,  p ). 

The  two  relations  (2.5a/b)  define  a  subset  of  possible  initial  points  (x ,  p )  of  £ .  For  a  closer  analysis  of  the 

"\ 

structure  of  this  set  we  introduce  the  orthogonal  projections 

P-  E  ->  L(  R",  R"),  P  (x ,  p)  R"  =  rge  DpF(x,  p )  ,  (x,p)  €  E  (2.6a) 

Q:  E  -»I(R\ R"),  Q(x,p)  =  In  -  P(x,p),  (x,p)  e  E.  (2.6b) 

These  projections  are  C '-functions  of  (x,  p).  In  order  to  see  this,  let  (x0,  po)  €  £  be  given  and  choose  an  ortho¬ 
normal  basis  [e^ . er )  of  rge  DpF(x0,  pQ).  Then  there  exist  w,  e  R"  such  that  et  -  DpF(x0,  p^)Wi  and  obvi¬ 
ously  for  (x,p)  near  (x0,  p0)  the  mappings  t).  .  r\t(x ,  p)  =  DpF(x , p)w(,  i  =  1 . r  are  of  class  Cl  and 

(Tli(x,  p),...,T|r(x,  p)}  are  linearly  independent.  Hence,  since  rank  DpF(x,p)  =  r,  these  vectors  form  a  basis  of 
rge  DpF(x,p).  By  applying  the  Gram-Schmidt  process  we  obtain  now  an  orthonormal  basis  of  the  same  space 
and,  because  the  process  involves  only  analytic  operations,  the  vectors  of  this  basis  arc  again  C1  functions  of 
(x.p).  But  then  the  same  holds  for  the  projection  />(x,  p)  as  the  sum  of  dyadic  products  of  the  basis  vectors 
whence  also  Q  =  /„  -  P  is  Cl. 

Evidently,  the  points  (x ,  p )  e  £  satisfying  the  necessary  conditions  (2.5a/b)  are  the  common  zeros  of  F 
and  the  C1  -mapping 

G  :  £  — »  R" ,  G(x,p)  =  />(x,p)£(x,p)  +  Q(x,p)D1£(x,p)p,(x,p)e  £.  (2.7) 

In  other  words, 

En  =  ((x.p)  e  £ ;  £ (x ,  p )  =  0,  G (x ,  p )  =  0)  (2.8) 

is  the  set  of  all  points  in  £  that  satisfy  the  necessary  conditions  of  Lemma  2.1.  For  later  use,  note  that  the  set 
( x(t ),  x' (0),  f  e  ]  lies  in  En  (and  not  only  in  £_1(0))  for  every  C2-soIution  of  (2.1). 

A  major  step  toward  transforming  (2.1)  locally  into  an  explicit  ordinary  differential  equation  will  be  pro¬ 
vided  by  the  following  result  about  the  relationship  between  the  solutions  of  (2.1)  and  those  of  the  equation 

C(x,x')  =  0.  (2.9) 

Below  and  in  other  instances,  we  shall  use  the  remark  that  G(x,  p)  =  0  amounts  to  the  fact  that  both  terms  on 
the  righthand  side  of  (2.7)  vanish. 


Lemma  2.2:  Any  C2-solution  (2.3)  of  (2.1)  solves  the  equation  (2.9),  and  conversely,  any  C2-solution  (2.3)  of 

(2.9)  for  which  F(x(t 0,  x' (t 0)  =  0  for  some  txs  J  is  a  C2-solution  of  (2.1). 

Proof:  By  Lemma  2.1  any  C2-soluiion  (2.3)  of  (2.1)  satisfies  F(x(t),  x'  (t))  =  0  and  (2.4)  for  all  t  €  7  whence, 
because  of  QDpF  =  0,  we  sec  that  Q(x,  x')DxF(x,  x')x'  =  0  on  all  of  7.  This  shows  that  x  is  a  C2-solution  of 

(2.9) . 

Conversely,  suppose  that  (2.3)  is  a  C2-solution  of  (2.9).  Then,  on  J  the  identity 
~ FQc,x ')  =  [P(x,x)  +  Q(x,x)][DxF(x,x')x'  +  DpF(x,x)x"] 

=  P(x,x')[DxF(x,x)x  +  DpF(x,x)x']=P(x,x')  4  F(.x,x')  , 

at 

holds  where  we  used  that  P  +  Q  =  /„,  QDpF  =  0,  and  Q  (x,  x')DxF(x,  x’)x!  =0.  By  differentiation  of  the 
identity  P  (x ,  x  )F  (x ,  x  )  =  0  (sec  (2.7)  and  (2.9))  it  follows  that 

P(x,xf)  4  F(x,x’)  =  -  4  P(x,x')  F(x,x').  Hence  the  function  J  ->R",^(i)  =  F(r(iU(i)) 
at  at 

satisfies  the  linear  system  4  ^  =  A(r)^,  where  A:  J  ->  L(R"),  A(r)  =  -  4  P(.x,x f)  ( t )  is  a  continuous 

at  at 

function  since  the  projection  P  is  C1  on  £  and  the  solution  x  is  C2  on  /.  Hence  the  standard  uniqueness  theory 
for  linear  systems  2)  together  with  the  condition  F(x(tx),  x'(tx))  =  0  implies  that  F (x  (t ),  /(r))  *  0  for  r  in  a 
neighborhood  of  ti  in  7.  Because  of  the  connectedness  of  7  it  follows  readily  that  this  local  result  holds  for  all 
t  €  7.  Thus  jc  is  a  solution  of  (2.1)  on  7  as  claimed. 

In  order  for  the  equation  (2.9)  to  induce  a  unique  vector  field  on  some  neighborhood  of  any  given  point 
(*o.  Po)  e  En,  we  need  to  guarantee  that  for  each  x  near  x0  there  exists  only  one  vector  p  near  p0  for  which 
G(x,  p)  =  0.  Obviously,  a  sufficient  condition  for  this  will  be  that  DpG(xo,  po)  is  an  isomorphism  on  R*. 
Accordingly,  we  define  the  set  of  admissible  initial  points  of  (1.1)  in  E  as 

Ea  =  [(x,  p)  e  En\  DpG(x,  p)  g  Isom( R")  }  .  (2.10) 


2)  We  emphasize  that  continuity  of  A  with  respect  to  t  is  sufficient  to  guarantee  uniqueness,  as  a  straightforward 
verification  confirms. 
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Clearly,  by  the  continuity  of  DPG,  the  set  EA  is  (relatively)  open  in  £w.  Moreover,  note  that  when 
(x0.  Po)  6  £.v  then,  by  (2.5b;,  there  exists  a  q0  e  R"  such  that 

DzF(x0,  po)Po  +  DpF(x0,  po)qo  =  0  .  (2.11) 

The  following  lemma  provides  a  characterization  of  EA  in  terms  of  F  and  its  derivatives.  For  a  somewhat 
related  condition  in  the  quasilinear  case  see  also  [2]. 

Lemma  2.3:  For  any  (jc0,  Po)  €  £.v  we  have  Dp  G  (x0,  p0)  e  Isom{  R")  if  and  only  if  for  some  q  0  which  satisfies 
(2.11)  the  following  condition  holds 

u  e  leer  DPF (x0,  p0 )  and  Dx2pF(x 0,  p0)(p0,  u)+DppF(x0,  p0)(q 0,  u)+DxF(x0,  p0)u  e  rge  DpF(x0,p o)  3) 
together  imply  that  u  =  0  .  (2.12) 

This  equivalence  does  not  depend  on  the  particular  choice  of  q0  satisfying  (2.11). 


Proof:  Evidently,  for  any  u  e  R  " ,  we  have 

DpG(x,  p)u  =  Kx{x,  p,  u)  +  K2(x,  p,  u) 


where 


Ki(x,p,  u)  =  Dp[P(x,p)F(x,p)]u  =  [DpP(x,p)u]F{x,  p) 

+  P{x,  p)DpF(x ,  p)u  -  [DpP(x,  p)u]F(x ,  p)  +  DpF(x,  p)u 

and 


K2(x,p,u)  *  Dp[Q(x,  p)DxF(x,  p)p]u. 

The  condition  (jc0.  Po)  e  implies  that  F(x0,  p0)  =  0  and  hence  by  (2.13)  that 
*■•  »  »  ■ 

K\(xo,P*  u)  =  DpF(x0,  Po)u. 


For  the  evaluation  of  K2  note  that  for  any  fixed  q  6  R"  and  by  definition  of  Q  it  follows  that 


Q(x ,  p)DxF(x , p)p  =Q(x,p)[DtF(x,p)p  +  DpF(x,p)q], 


and  therefore  that 

3)  Here  DxpF(x0,  Po)(po>  a)  mean*  D2F (x0,  Po)((Po.O),(0,M ))  »nd  hence />0  »nd  u  do  no*  ptsy  iyran»e«ric 
rotes. 


(2.13) 


(2.14) 


(2-15) 


(2.16) 
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K2(x,p,  u)  =  {Dp Q(x,  p)u ){DX F (x,  p )p  +  DpF(x,  p)q) 

+  Q(X.  p)[Dz2pF(x,  p)(p,  u )  +  DxF(x,p)u  +  DppF (x ,  p)(q ,  u) ]. 

Now  with  any  point  qQ  e  R'1  for  which  (2.1 1)  holds  we  find  that 

i<2(xo,Po,u)  =  Q(x0,  Po)[DxzpF {xQ,  p0)(po,  u)+D2pF (x0,  p0)(qQ,  u}+DxF (x0,  p0)u]. 

Together,  (2.15)  and  (2.16)  show  that  Ki(x0,  p0,  u)  and  K2(x0,  p0,  u)  are  the  components  of  the  vector 
DpG(x0,po)u  along  rge  DpF (x0,  p0)  and  its  orthogonal  complement,  respectively.  This  implies  that 
DpG(x0,  p0)u  =  0  if  and  only  if 

K^xq.pq,  u)  =  0  and  K2(x0,  p0,  u)  =  0.  (2.17) 

Thus  DPG  (x0,  pQ)  is  an  isomorphism  exactly  if  (2.17)  holds  only  for  u  =  0  which  by  (2.15)  and  (2.16)  is 
equivalent  with  (2.12).  The  last  part  of  the  lemma  now  is  a  direct  consequence  of  the  fact  that  the  invertibility 
of  DpG(x0 ,  p 0)  is  independent  of  the  choice  of  q0. 

As  an  immediate  corollary  of  Lemma  2.3  we  obtain  from  (2.10)  and  the  implicit  function  theorem  the  fol¬ 
lowing  result: 

Lemma  2.4:  For  any  ( x0 ,  p0)  e  EA  there  exists  an  open  neighborhood  UQ  x  V0  c  E  and  a  unique  Catnapping 
<I> :  U0  -»  y0,  ^(jco)  =  Po<  such  that  (x ,  p )  6  U0xV0  and  G (x ,  p)  =  0  if  and  only  if  p  =  <&(x). 

Lemma  2.4  shows  that  for  any  initial  point  (x0,  Po) 6  Ea  there  exists  an  open  neighborhood 
E0  =  UqX  V0  in  E  where  the  system  (2.9)  can  be  written  in  the  explicit  form  x'  =  <P(x).  Hence  the  standard 
theory  for  initial  value  problems  ensures  that,  modulo  translations  in  time,  this  explicit  system  has  a  unique  solu¬ 
tion  x  in  C0  through  any  given  point  of  that  set,  and  clearly  this  solution  is  of  class  C2  since  O  is  C1.  It  follows 
that  for  any  given  point  of  £o  the  system  (2.9)  has  a  unique  C2  solution  x  such  that  (x,  x')  passes  through  that 
point  of  £0.  Now  Lemma  2.2  asserts  that  such  a  solution  is  a  solution  of  (2.1)  if  and  only  if  (x.x')  passes 
through  some  point  of  £ Thus  we  conclude  that  (2.1)  has  a  unique  C2  solution  for  which  (x.x')  passes 
through  any  given  point  of  EN  ^  £0.  As  noted  earlier,  EA  is  open  in  EN  and,  clearly,  for  sufficiently  small  £ 0 
we  have  EN  pi  £0  =  EA  (~\E0.  Thus,  in  particular,  the  result  applies  to  the  given  point  (x0,  po)  e  EA. 
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We  summarize  this  conclusion  in  the  following  form: 

Theorem  2.1:  Suppose  that  for  the  problem  (2.1)  the  conditions  (2.2a/b)  arc  valid  and  choose  any  r0  e  R .  Then, 
(*o.  P o)  e  £,v  is  necessary  for  the  existence  of  a  C2-solution  of  the  initial  value  problem 

F  (x ,  x  )  =  0,  x(t0)  =  jc0.  x'  (to)  =  p0  .  (2.19) 

Moreover,  if  Cc0,  p0)  €  EA  then  there  exists  a  unique  C2-solution  of  (2.19). 

Note  that  the  theorem  does  remain  valid  when  r  =  n  in  (2.2b).  In  that  case  we  have  rge  DpF(x,p)  -  R* 
for  all  (x ,  p )  e  E  and  hence  the  necessary  condition  (x0,  p q)  €  EN  simply  is  F(x0,  p0)  =  0  while  the  sufficient 
condition  (x0,  p 0)  e  Ea  reduces  to  F(x0,  p0)  =  0  and  ker  DpF(xa  ,  p0)  =  (0).  Thus  we  recover  here  the  usual 
situation  when  the  implicit  function  theorem  provides  that  (2.1)  can  be  written  locally  as  an  explicit  ODE.  In 
this  case,  it  is  of  course  sufficient  that  F  be  of  class  C1  on  £  since  the  second  derivatives  of  F  are  no  longer 
involved  in  the  definition  of  EA . 

Similarly,  the  extreme  case  r  =  0  in  (2.2b)  is  trivial  since  then  DpF(x,  p)  »  0  for  (x,  p)  e  E ;  that  is,  F 
is  independent  of  p.  Here  we  have  (x0,  pQ)  e  EN  if  and  only  if  F(x o)  =  0  and  DxF(x0)po  =  0  while 
(xo,  po)  €  Ea  under  the  additional  requirement  that  DXF  is  at  x0  an  isomorphism  of  R*  to  itself.  But  then  the 
equation  F  =  0  has  x0  as  isolated  solution  and  the  unique  solution  of  (2. la/b)  is  x(r)  =  xo.  This  is  consistent 
with  the  remark  that  if  DxF(x0)  is  an  isomorphism,  then  DxF(x0)p0  =  0  only  if  p0  =  0. 

We  end  this  section  by  considering  the  conditions  (x0,  Po)  e  Es  and  (x0,  p0)  e  EA  of  Theorem  2.1  for 
the  general  nonautonomous  case  (1.1);  that  is, 

F(t,  x,  xf)  =  0,  x(t0)  =  xo,x'(ro)  =  po,  F(fo,Xo,Po)  =  0.  (2.20) 

We  use  the  standard  approach  to  make  this  problem  autonomous  and  hence  introduce  the  mapping 

■  ’ 

//:  (Rx  R")x(Rx  R")  — »  R*+l  ,  // ((r ,  x),  (t,  p))  =  Xp[t)Xt{})  (2.21) 

and  the  corresponding  initial  point  ((r0»  x0),  (t0,  p0)),  t0  =  1.  Then,  under  the  required  smoothness  assumptions 


the  condition  (2.5b)  for  H,  written  in  terms  of  F,  assumes  the  form 

D,F°  +  DxF°p0  €  rge  DpF° 


(2.22) 
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while  the  condition  (2.12)  becomes 

for  q 0  e  R"  such  that  DlF°+D1F0p0+DpF0q0  =  0, 
u  €  ker  DpF°  and 

D,^F°u+  D^F°(po,  u)+DppF°(q0,  u)+DxF°u  e  rge  DpF°  (2.23) 

together  imply  that  u  =  0  . 

In  (2.22)  and  (2.23)  the  superscript  0  indicates  that  the  particular  function  is  to  be  evaluated  at  (r0,  xq,  po ).  Note 
also  that  in  (2.22)  and  (2.23)  we  explicitly  used  t0  =  1  and  that  the  condition  (2.23)  is  independent  of  the 
specific  choice  of  q0. 


3.  Global  Behavior  of  the  Solutions 

In  order  to  determine  the  global  behavior  of  the  local  solutions  of  the  initial  value  problem  (2.19)  suppose 
again  that  the  conditions  (2.2a/b)  hold.  Then  Theorem  2.1  guarantees  the  existence  of  a  unique  C2-solution  (2.3) 
x  for  which  (x ,  x' )  passes  through  any  given  point  of  EA ,  For  any  such  solution  and  any  ye  )  we  introduce  the 
sets  T^-  [p  g  R“;  p  =  x'  (t),  y  St  <  b }  and  Tr  =  [p  e  R";  p  =  x'(r).  a  <  t  gy).  Then  the  following 
extendability  result  holds: 

Theorem  3.1:  Assume  that  for  the  problem  (2.1)  the  conditions  (2.2a/b)  hold  and  that  EA  =  EN.  Then  the  follow¬ 
ing  statements  are  valid: 

(i)  If  the  set  T y*  [or  Ty.]  is  bounded  for  some  ye/,  then  lim,_^_  x(t)  =  xb  [or  lim,_*,+  x(r)  =  x.]  exists. 

(ii)  If  lim,_*_  x(r)  =  xb  [or  lim,_KI+  x(t)  -  xj  exists  and  for  some  sequence  [rt )  e  J  with  limt_*.ft  =  b  [or 
limt^r*  =  a]  the  sequence  [x'(r*))  has  an  accumulation  point  p*  for  which  (xb,p*)eE  [or 
(xa,  p*)  e  E]  then  lim,_*_x'(f)  =  p*  [or  lim^^  x' (t)  =  p*)  and  for  b  <  «o  [or  a  >  the  solution  x 
can  be  continued  to  the  right  of  b  [or  to  the  left  of  a]: 

Proof:  We  present  the  proof  only  for  the  right  endpoint  b,  for  the  other  one  it  proceeds  analogously.  Suppose 
first  that  7V  is  bounded  for  some  y  g  /  and  hence  that  llx'  (r)ll  <  M  <  <*>  for  t  e  /T  =  [y,  b).  Since 

I 

llx(r)  -  x(j)ll  S  I  J  llx'(t)lldxl  £  Af  Ir-s  I  ,  for  aU  s,  t  e  /T 


-  10  - 


we  see  that  for  any  sequence  {4}  c  with  limk_rsr-^b  the  sequence  [x(tk)}  is  a  Cauchy  sequence,  and 
moreover  that  its  limit  is  independent  of  Lhe  particular  choice  of  {/*}.  This  proves  the  existence  of  the  limit 
point  xb . 

Now  assume  that  lim^j,..  x(r)  =  xb  exists  and  that  for  some  sequence  [tk]  <i  J  with  limt_*_/t  =  b  the 
sequence  {x'(rt)}  has  an  accumulation  point  p*  for  which  (xb,p*)e  £.  Then  there  is  a  sub-sequence 
[sk )  c  {t*}  for  which  lim*.,^*  =  b  as  well  as  limt__r'(jt)  =  p * ,  and,  because  of  F(x(sk),  x'  (sk))  =  0  and 
(xb ,  p* )  e  £ ,  we  also  have  F(xb,  p*)  =  0.  As  in  the  previous  section,  for  all  (x ,  p )  near  (xb,p*)  let  Q (r ,  p ) 
be  the  orthogonal  projection  onto  [rge  DpF(x,p)}  .  Recall  that  Q  is  continuous  (even  C').  Hence,  from 
DxF{x(sk ),  x  (sk))x  (sk)  e  rge  DpF (x (sk),  x  (sk ))  (see  Lemma  2.1)  it  follows  that  Q(x(sk),  xf  ( sk))DMF(x(sk ), 
x’  (sk))x!  (r*)  =  0  for  all  k  and  therefore,  in  the  limit,  that  Q  (xb ,  p*  )DXF (xb ,  p*  )p*  =  0.  In  other  words,  p*  is  a 
solution  of  the  system 

F(xb,p)  =  0  ,  DxF(xb,p)p  e  rge  DpF  (xb,  p)  .  (3  2) 

Since  (xb,p*)  e  E  it  follows  that  {xb,p*)  e  EN  =  EA  and  thus,  by  Lemma  2.3,  that  the  solution  p  =  p*  of 
(3.2)  is  isolated.  In  other  words,  there  exists  a  ball  BjcR"  centered  at  p*  with  radius  8>0  such  that 
(xb,  p)  e  £ ,  for  all  p  e  Bb,  and  that  fl5  contains  no  solution  p  of  (3.2)  other  than  p* . 

Consider  now  any  sequence  (x* }  e  J  with  lim*_^t*  =  b .  We  show  first  that  x' (t*)  must  have  an  accu¬ 
mulation  point.  Suppose  that  this  is  not  true,  so  that  limJk_^.Ilx,(xi)l!  =  «=.  Because  of  limt_+-x'(jt)  =  p*  it  there¬ 
fore  follows  that  for  some  sufficiently  large  kQ 

IUt'  (Tt)  -  *'  (j*)II  >  5/2,  for  k  £  k0.  (3.3) 

By  continuity  and  the  intermediate  value  theorem  there  exists  in  each  interval  [inf{xk,  sk),  sup (x* ,  s*)]  a  value 
p*  such  that  Hx'(p*)  -  x  (sk )ll  =  5/2.  Then,  x  (p* )  is  bounded  and  upon  extracting  a  subsequence,  we  obtain  an 
accumulation  point  p p  of  x'(p4),  which,  of  course,  satisfies  llpp  -  p*  il  =  6/2.  But  this  means  that  pp  e  B j  and 
*  p*  which  is  a  contradiction  because  (xb ,  p p)  e  £,  limt^p*  =6  and  hence  pp  is  a  solution  of  (3.2). 
Therefore  x  (r* )  must  have  an  accumulation  point.  Suppose  there  exists  such  an  accumulation  point  p**  which 
is  distinct  from  p* .  If  ( xb ,  p**)  e  £  then,  as  before,  p **  must  solve  (3.2)  and  hence  cannot  be  in  Bj.  On  the 
other  hand,  if  ( xb ,  p** )  docs  not  belong  to  £  then  necessarily,  p**  cannot  be  in  £5  either.  There  is  some  sub- 
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sequence,  to  be  denoted  again  by  {t*  } .  which  converges  to  b  and  for  which  lim*_*.  x'Ct*)  =  p** .  As  p**  is 
not  in  B&,  the  relation  (3.3)  must  be  valid  for  sufficiently  large  k0  which,  of  course,  leads  once  more  to  a  con¬ 
tradiction.  Thus,  altogether,  we  have  found  that  p*  is  the  unique  accumulation  point  of  the  sequence  {x'Cr*)) 
for  any  sequence  {t* }  e  J  such  that  lim*_,.  zt  =  b .  This  proves  that  lim,.,*,  x'(f)  =  p* . 

Since  p*  solves  (3.2)  we  have  ( xb ,  p*)  e  EN  and,  by  hypothesis,  EN  =  EA.  Hence,  for  b  <  <»  we  can 
apply  Theorem  2.1  at  (x0,  pQ)  =  (xb ,  p* )  with  t0  =  b.  Thus  there  exists  some  neighborhood  £0  c  E  of  (xb ,  p* ) 
and  some  open  interval  /  containing  b  where  the  C2-solution  y:  /  -»  R"  y(b)  =  xb,  y’(b)  =  p* ,  of  (2.1a)  is 
the  only  solution  of  an  explicit  system  y  =  d>(y),  x(b)  =  xh,  with  some  C '-function  <t>  on  E0  for  which 
<&(*»)  =  p* .  By  the  standard  uniqueness  theorem  it  follows  that  any  "one-sided"  C '-solution 
{b  -  e,  b]  -»  R"  with  (q(t),  q'(O)  e  E0  ,  b  -  e  <  t  <  b,  of  this  explicit  problem  necessarily  has  to  agree 
with  the  unique  C2-solution  y  of  the  problem  on  their  common  interval  of  definition.  But  because  of 
lim,_^_  x(t)  =  xb,  and  lim,_^_  x'(t)  =  p*  our  given  solution  x  is  such  a  one-sided  C'-solution  and  hence 
agrees  with  y  on  their  common  domain.  This  shows  that  the  original  solution  indeed  can  be  continued  beyond 
the  right  endpoint  b  of  J  and  the  proof  is  complete. 

The  result  implies  that  any  local  C2-solution  of  (2.1a)  can  be  extended  to  some  open  interval 
J  =  (a* ,  b* ),  b*  <,  a*  >  which  is  maximal  under  set  inclusion. 

Now  consider  Theorem  3.1  in  the  case  £=R'lxR’,.If  limx_fc_x(r)  =  xb  exists,  then  any  accumulation 
point  p  of  x(t)  as  t  -♦/)-  in  R"  necessarily  satisfies  (xb,  p)  e  E  and  we  have  lim,_^_  x'(f)  =  p.  Clearly,  if 
is  bounded  then  there  must  be  such  an  accumulation  point;  in  other  words,  for  bounded  both  limits 
lin**-**--^*)  and  limx^*_x' ( t )  always  exist  and  the  solution  can  be  extended.  Suppose  now  that  /»=/>•<« o 
which  implies  that  the  solution  cannot  be  extended.  Then  7^  must  be  unbounded.  Moreover,  if 
limx_^_  x(t)  =  xb  exists  then  x' (t)  cannot  have  an  accumulation  point  a s  t  -*  b-  which  means  that 
limx^*_  IU'(t)ll  =  °°.  The  analogous  result  holds  at  a*  and  altogether  we  have  the  following  corollary  of 
Theorem  3.1. 

Theorem  3.2:  Assume  that  the  mapping  F  of  (2.1a)  satisfies  the  conditions  (2.2a/b)  and  (3.1)  on  £  =R*xR" 
and  let  x:  /-*  R*.  (x(0,  ^ (f))  6  £  for  t  e  J ,  be  any  C2-solution  of  (2.1a)  where  /  »  (a*,  b*)  is  a  maximal 
interval.  If  b*  <  *>  then  the  set  is  unbounded  for  some,  and  hence  all,  yeJ,  and,  in  particular,  if 
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x(f )  exists  then  lim,^..  ILc'(t)li  =  Correspondingly,  for  a*  >  -«o  the  set  Ty.  is  unbounded  and,  in 
particular,  if  lim,_«.+  x(t)  exists  then  lim,_0.+  lbc'(r)ll  = 


4.  The  Structure  of  the  Set  of  Admissible  Initial  Points 

In  this  section  we  analyze  the  structure  of  the  set  EA  of  admissible  initial  points  of  (2.1)  in  £  as  well  as 
that  of  its  projection  onto  the  x-spacc.  For  this,  we  assume  that,  in  addition  to  the  properties  (2.2a/b)  the  map¬ 
ping  F  satisfies 

rank  DF(x,  p)  =  n,  for  all  (jc,  p)  e  £;  (4.1) 

that  is,  that  the  equations  (2.1)  arc  independent. 

We  prove  first  the  following  preliminary  lemma: 

Lemma  4.1:  Suppose  that  F  satisfies  the  conditions  (2.2a/b)  and  (4.1).  Then  the  mapping  QF :  E  -»  R"  has  at 
any  point  (x,  p)  e  £_l(0)  the  partial  derivatives 

Dx(QF)(x,p)  =  Q{x,p)DxF(x,p),  Dp(QF)(x,p)  =  0.  (4.2) 

Moreover,  the  linear  map  Q(x,  p)DxF(x,  p)  €  L(R")  has  constant  rank  n-r  on  all  of  F-1(0). 

Proof:  Let  (x,  p)  e  £-1(0).  For  any  u  e  R"  we  have 

Dx[Q{x,p)F(x,p)]u  =  [DzQ(x,p)u]F(x,p)+Q(x,p)DxF(x,p)u 

which  implies  the  first  part  of  (4.2)  while,  because  of  QDPF  ■  0,  the  second  part  is  a  consequence  of 

Dp  IQ  (x ,  p  )F  (x ,  p )]  u  =  [Dp  Q  (x ,  p  )u  ]F  (x ,  p  )+Q  (x ,  p  )Dp  F  (x ,  p  )u. 

By  (4.2)  we  see  that  Q(f.p)DxF(x,p)  =  Q(x,p)DF(x,p)  and,  because  DF(x,p)  has  full  rank,  that 
rank  Q(x,  p )DF (x ,  p )  =  rank  Q(x ,  p)  =  n-r .  This  proves  the  second  part  of  the  assertion. 

With  this  we  obtain  now  the  following  result  about  the  structure  of  EA : 

Theorem  4.1:  Suppose  that  F  satisfies  the  conditions  (2.2a/b)  and  (4.1).  Then  the  set  EA  cE  of  admissible  ini¬ 
tial  points  of  (2.1)  is  an  r -dimensional  Cl -sub-manifold  of  R"  x  R*. 
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Proof:  Let  (x0,  Po)  e  EA ,  so  that  F(x0,  p0)  =  G(x0,  p0)  =  0  and  DpG(x0,p Q)  e  Isom( R").  Since  the  condition 
DpG(x,p)  e  Isom (R")  is  fulfilled  by  all  points  (x ,  p )  e  E  near  (jc0,  p0),  the  points  of  EA  in  the  vicinity  of 
C*o»  Po)  characterized  by  the  sole  conditions  F(x,  p)  =  G(x ,  p)  =  0.  Actually,  this  system  may  be  replaced 

by 


G*(x,p)  =  (Q  (x0,  Po)F(x,p)  ,  G  (x.  p))  =  0  . 


Indeed,  it  is  obvious  that  G*(x,p)  =  0  whenever  F(x,p)  =  G(x,p)  =  0.  Conversely,  suppose  that 
G*(x,p)  =  0,  so  that  Q(xq,  p0)F(x,  p)  =  0,  G(x,p)  =  0.  The  latter  relation  implies  in  particular  that 
P(x,p)F(x,p)  =  0,  which,  together  with  Q(x0,  po)F(x,  p)  =  0  yields  F(x,p)  =  0  because 
rg<  Q(x,p)  =  [rge  DpF {x,  p)Y  and  rge  P  (x0.  p0)  ~  rge  DpF(x0,  p0)  remain  complementary  for  (x,  p)  close 
enough  to  (xo.  Po)  (by  constancy  of  rank  DPF  and  continuity  of  P  and  Q  on  £;  see  Section  2). 

Identifying  rge  Q(x0,  p0)  =  R’,"r,  we  see  that  the  mapping  G*  maps  a  neighborhood  of  (xo.  Po)  in 
R"x  R"  into  the  fixed  space  R"_rx  R",  and  with  F(x0,p o)  =  0  and  QDpF  =  0,  that 


DG*(x0,p  q)  = 


Q(xo,  Po)F>*F(x0,  po) 
PxG (x0,  po) 


0 

DpG(xq,  Pa) ' 


(4.3) 


Recalling  that  DpG(x0,  p0)  e  Isom( R")  and  from  Lemma  4.1  with  (x,p)  =  (x0, Pa),  we  conclude  that 
DpG*(x0,  po)  maps  onto  R"'rxR",  The  implicit  function  theorem  now  ensures  that  (G* )-I(0),  and  hence  EA 
is  a  r-dimensional  Cl  sub-manifold  of  R"x  R"  in  the  vicinity  of  (xo. po).  This  completes  the  proof  since 
(xo,  Po)  was  an  arbitrary  point  of  EA . 

Let  n  :  R"x  R*  -»R"  be  the  projection  onto  the  first  factor  and  let  (x0,  po)  e  EA.  It  follows  from  the 
proof  of  Theorem  4.1  that  the  tangent  space  T{lyP^EA  may  be  identified  with  the  null-space  of  the  mapping 

DG *  (x0.  Pa)  in  (4.3).  Thus,  for  given  (u ,  q )  6  T p^EA ,  one  has  n(u  ,q)  =  u  and  hence  n(u ,  q )  =  0  if  and 
only  if  («,<?)  =  (0,  0).  This  means  that  the  restriction  of  n  to  the  manifold  EA  is  an  immersion  at  (xo,  Po).  It 
follows  that  there  is  an  open  neighborhood  N0  of  (x0,  Po)  in  EA  such  that  Af0  *  n(No)  is  an  r-dimensional  Cl 
sub-manifold  of  Rn  and  U.  :  N0~*  Mo  *s  a  C 1  -diffeomorphism.  We  are  now  in  a  position  to  prove  that, 
locally,  the  DAE  (2.1)  is  equivalent  to  an  explicit  ODE  on  the  manifold  M0 


Theorem  4.2:  There  is  a  neighborhood  E0  =  U0  x  V'q  c  R*x  R"  of  (x0,  Po)  e  EA  and  a  Cl  vector,  field 
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O  :  U0  Mq  ->  Rn ,  <t>(x0)  =  Pa.  with  <b(x)  e  TXM0,  for  x  e  U0  C\  M0,  such  that  x  :  J  -*  R*  is  a  C2- 
solution  of  E(x.x')  =  0  satisfying  (x(r),  x'(t ))  e  E0,  for  teJ,  if  and  only  if  x(t)eU0r^Mo- 
x'(t)  =  <&(x(t )),  for  all  t  s  / . 

fYoo/:  Lemma  2.4  has  already  established  existence  of  a  neighborhood  E0  =  UqxVq  and  C'-mapping 
d> :  U0  Rn,  <b(x0)  =  p o,  such  that  (x,  p)  e  E0  and  G(x,  p)  =  0  if  and  only  if  x  e  U0  and  p  =  d>(x).  More¬ 
over,  it  was  also  shown  that  x:  J  R"  is  a  C2  solution  of  (2.1)  satisfying  ( x(t ),  x'(t))  6  £0  if  and  only  if 
x(r)  e  t/0.  x'(r)  =  d>(jt(r)),  for  all  r  €  7,  and  (x,  x')  passes  through  one  point  of  E*  (so  that,  in  fact,  (x,x') 
lies  entirely  in  £\v). 

Clearly,  the  neighborhood  E0  =  U0x  V0  can  be  shrunk  to  arbitrarily  small  size  with  no  prejudice  to  the 
mentioned  properties.  One  may  then  assume  that  N0  =  E0  fy}  EA  =  £0  E"  in  the  discussion  preceding  the 
theorem.  In  this  case,  M 0  =  n(JV 0)  c  U0  so  that  U0fy^M0  =  M0. 

Since  the  curve  (x,  x')  lies  in  EN  whenever  x  is  a  C2-solution  of  E(x,  x)  =  0,  it  follows  that  x(f)  6  M0 
as  soon  as  (x(r),  x'(r))  e  £0.  To  prove  the  theorem,  it  suffices  to  show  that 

(a)  x  e  M0  implies  (x,  d>(x))  e  N0, 

(b)  x  e  M0  implies  d»(x)  €  TXM 0  c  R". 

If  (b)  is  assumed  to  hold  then,  (a)  is  needed  to  show  that  for  any  (automatically  C2)  solution  x:  J  -*  M0 
to  x'  =  <D(x)  the  curve  (x,  x')  lies  in  N 0  c  EN  and  hence  is  a  C2-solution  of  (2.1).  But  (a)  holds  since  every 
x  e  Mq  has  the  form  x  =  n(x,  p)  with  (x,  p)  e  N0  c  £0,  whence  p  =  <t>(x). 

To  prove  (b),  we  shall  use  a  characterization  of  the  points  of  N0  that  is  slightly  different  from  that 
involved  in  the  proof  of  Theorem  4.1:  Note  that  (x ,  p )  e  N0  if  and  only  if  Q(x,p)E(x,p)  =  0  and 
C(x,p)  =  0  since  £(x,p)£(x,p)  =  0  is  already  ensured  by  the  second  relation  and  since 
DpG(x,  p)  €  /rom(R")  is  guaranteed  by  the  hypothesis  E0  EN  =  E 0  E*  (=#o)-  Since  (x.p)  e  £<>,  we 
see  that  G(x,  p)  =  0  is  equivalent  with  x  e  C/0  and  p  =  4>(x)  and  hence  (x.p)  6  Af0  if  and  only  if  x  e  U0  and 
H  (x )  =  0  where  we  have  set 

//(x)  ■  Q(x,  <t>(x))E(x,  <t>CO)  ,x  e  f/0  • 

Using  (a),  one  finds  at  once  that  Ma  =  //_1( 0).  In  turn,  this  shows  that  for  x  e  A/0.  the  space  TmMq  identifies 
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with  ker  DH(x).  Because  of  //(x)  =  0  we  have  (x ,  <t>(x))  e  N0  c  EN  c  E_1(0)  and,  with  the  help  of  Lemma 
4.1,  an  elementary  calculation  yields 

DH(x)u  =  Q(x,  $>(x))DxF(x,  <b(x))u  , 

for  u  e  R\  Letting  u  =  <t>(x),  we  find  that  <I>(x)  e  ker  DH(x)  =  TXM0  because 
Q(x,  <t>(x))DzF(x,  <t>(x))d»(*)  =  0  follows  from  G(x,  <t>(x))  =  0.  This  completes  the  proof  of  Theorem  4.2. 

Remark:  In  the  proof  of  Theorem  4.1,  the  projection  Q(x0,  p0)  cannot  be  replaced  by  Q(x,p),  for  rge  QQc,p) 
varies  with  (x,  p)  and  the  resulting  mapping  G*  must  be  viewed  as  taking  its  values  in  R*x  R*  instead  of 
R*~,x  R",  and  the  implicit  function  theorem  cannot  be  used.  Neither  can  the  rank  theorem  be  applied  since 
there  is  no  guarantee  that  rank  DG*(x,  p)  =  2n-r  for  (x,p)  near  (x0,p q)  but  not  on  EA.  On  the  other  hand, 
Q(x0,  po)  cannot  be  substituted  for  Q(x,p)  in  the  proof  of  Theorem  4.2,  for  it  would  become  impossible  to 
take  advantage  of  Lemma  4.1  with  (x,  p)  *  (x0,  po). 

As  in  [7]  the  local  result  in  Theorem  4.2  may  be  globalized  to  some  extent  by  applying  the  theory  of  cov¬ 
ering  spaces.  We  sketch  only  briefly  the  general  approach.  Clearly,  the  local  result  shows  that  the  restriction 
nA  *  n  \Ea  is  a  local  homcomorphism  between  EA  and  nA£A .  Let  £A  be  some  non-empty,  arc-connected  sub¬ 
set  of  Ea  for  which  (EA,  nA),  nA  =  n l£A,  is  a  covering  space  of  FI£A.  In  other  words,  each  point  x  e  n£A  is 
assumed  to  have  an  open,  arc-connected  neighborhood  U  such  that  each  arc-component  of  (IIA)_1(C/)  is  not 
empty  and  is  mapped  topologically  onto  U  by  nA.  Often  it  turns  out  that  £A  =  £A  can  be  used.  This  is  certainly 
the  case  when,  for  Axed  x  e  n£A  there  are  only  finitely  many  p  such  that  (x,  p)  €  £A.  In  general,  it  is  always 
possible  to  choose  £A  as  the  closure  of  a  non-empty,  pre-compact,  (relatively)  open,  and  arc-connected  submani¬ 
fold  of  £a  . 

For  any  given  (x0,  p  o)  e  £A  let  now  W  be  a  non-empty,  (relatively)  open,  simply  connected  subset  of 
n£A  that  contains  xq.  For  any  x  e  M  choose  a  path  \  :  J  -*  M  which  connects  x0  with  x.  Then  there  exists  a 
unique  lifting  J  -»  £A  with  initial  point  (x0,  p0)  for  which  nA£*  =  £.  This  lifted  path  has  a  unique  endpoint 
(x,p)  in  £a  since  all  paths  in  M  between  x0  and  x  are  homotopic.  Since  x  was  arbitrary  in  M  our  above  local 
result  can  now  be  used  to  prove  that  M  indeed  is  an  r-submanifold  of  R*  and  that  the  DAE  (2.1)  induces  a 
tangential  vector  field  on  M  for  which  all  integral  curves  in  M  are  solutions  of  (2.1). 


i 


It  is  well-known  (see  c.g.  [1])  that  (5.1)  is  uniquely  solvable  for  compatible  initial  data  if  and  only  if  the 
matrix  pencil  (A ,  B )  is  regular;  that  is,  exactly  if  there  is  some  X.  e  R  such  that  B  +  XA  is  invertible.  A  central 
concept  in  the  solvability  theory  of  (5.1)  is  the  index  of  regular  pencils.  For  regular  ( A,B )  choose  any  X  for 
which  B  +XA  is  invertible.  Then,  the  index  is  the  smallest  integer  k  (£  n)  such  that 

ker  [(B  +  U  ylA  P1  =  ker  [{B  +  XA  )~'A  ]k.  (55) 

It  can  be  shown  that  k  is  independent  of  the  choice  of  X  (see  [4,  App.  A]),  and  it  is  also  readily  checked  that 
k  =  0  if  and  only  of  A  is  invertible. 

In  order  to  relate  our  theory  to  this  index-concept,  let  F  e  L(R")  be  the  orthogonal  projection  onto 
rge  A  and  Q  =  /„  -  P .  As  in  Section  2  our  first  step  is  to  differentiate  the  DAE  (5.1)  and  then  to  multiply  thr 
resulting  equation  Ax"  +  Bx'  -  f  by  Q  in  order  to  remove  again  the  second  derivative  of  x.  Together  with  the 
projection  of  the  original  equation  onto  rge  A  this  produces  the  system 

FAX'  +  PBx  =  Pf ,  QBx'  =  Qf  .  (5.6) 

Since  P  and  Q  map  onto  complementary  spaces  the  two  equations  (5.6)  can  be  added  which  results  in  die 


reduced  DAE 
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A  {x  +  B  [X  =  /  [  (5.7a) 

where,  because  of  FA  =  A , 

Ai=PA  +QB  =A  +  QB,  B  ,  =  PB  ,  / 1  =  Pf  +  Qf'  .  (5.7b) 

In  Section  2,  (5.7a/b)  corresponds  to  (the  non-autonomous  version  of)  (2.9)  and  thus  the  results  of  that  section 
can  be  applied.  We  concentrate  here  only  on  the  effect  of  the  sufficient  condition  (5.4).  Even  without  recourse  to 
the  earlier  theory,  it  is  readily  checked  that  (5.4)  is  equivalent  with  A,  e  Isom( R")  and,  hence,  that  when  (5.4) 
holds  then  (5.7a/b)  can  be  transformed  into  an  explicit  ODE. 

Suppose  therefore  that  A  {  is  singular.  Then  we  may  apply  the  same  procedure  repeatedly  as  many  times  as 
necessary,  to  obtain  a  sequence  of  DAEs  of  the  form 

Ajx’  +  BjX  =  f  j  (5.8a) 

wher eAj.Bj.fj  are  specified  recursively  by  A0  =  A,B0  =  B,f0  =  f  and 

A,+i  =  P jAj  +  Q,Bj,  Bj+i  =  PjBj,  /,+1  =  P,fj  +  Qjf'j ,  (5.8b) 

and  Pj  is  the  orthogonal  projection  onto  rge  Aj  and  Qj  =  /,  -  Pr  The  process  stops  with  the  smallest  integer  k 
such  that  Ak+l  is  invertible;  that  is  when  the  sufficient  condition  (5.4)  holds  for  A*  and  Bk.  Explicitly,  this  condi¬ 
tion  has  the  form 

Aku0  =  0  and  Bku0e  rge  Ak  imply  u0  =  0.  (5.9) 

In  terms  of  the  original  matrices  A  and  B  of  (5.1)  the  condition  (5.9)  turns  out  to  be  equivalent  with  the 
condition 

Auq  =  0,  Buj  =  Aiij+U  j  =0,1 k,  implies  u0  =  0.  (5.10) 

The  proof  will  follow  by  repeated  application  of  the  following  result  where,  for  ease  of  notation,  the  matrices 
A ,  B  now  stand  for  any  A, ,  BL : 

Lemma  5.1:  Let  A,  B  be  any  nxn  matrices  and,  with  the  orthogonal  projections  P  and  Q  ■  Ia-P  onto  rge  A 
and  ( rge  A  )J,  respectively,  set  A  i  =  PA  +  QB  =  A  +  QB ,  B\  =  PB .  Then,  for  any  k  £  1  the  equations 


Au0-0,  Buj  =  Auj+i,  j=Q  »•••» 


(5.11) 
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have  a  solution  uq,  1  if  and  only  if  u0 . uk  solve  the  equations 

A(u0  =  0,  B iiij  =  A  [UJ+l,  j=0 . *-l.  (5.12) 

In  particular,  if  for  all  solutions  of  (5.11)  we  have  u0  =  0,  then  the  same  must  hold  for  all  solutions  of  (5.12). 

Proof:  Suppose  that  u0,  ...,w*+i  solve  (5.11).  Since  QA  =  0,  multiplication  of  Buj  =  AuJk. x  by  Q  shows  that 
QBuj  =  0  for  j  =  0 Hence,  for  j  =  0,...^t-l  we  obtain 

BiUj  =  PBuj  =  PAuj+x  =  PAuJ+i  +  QBUj+x  =  A^ 

and,  since  Au0  =  0, 

Axu0  -  PAu0  +  QBu0  =  PAuo  =  0, 

so  that  Uq . uk  solve  (5.12). 

Conversely,  assume  that  the  vectors  u0,  ...,uk,  solve  (5.12).  In  terms  of  A,  B  the  equations  (5.12)  assume 
the  form 

Aa0  +  QBu0  =  0  (5.13a) 

PBuj  =Auj+1  +  QBuj+x,  j  =  0,...,  fc-1  (5.13b) 

where  we  used  that  A  t  =  A  +  QB .  Since  Q  maps  onto  a  complement  of  rge  A ,  both  terms  on  the  left  of  (5.13a) 
have  to  be  zero;  that  is,  we  have  Au0  =  0,  and  QBu0  =  0.  By  multiplication  of  (5.13b)  with  Q  it  follows  that 
QBuj+x  =  0  and  hence  altogether  that 

Buj+x  =  PBuj+i,  j  =  0,...,  *-1.  (5.14) 

Now,  by  multiplying  (5.13b)  with  P  and  using  (5.14),  we  obtain  Buf  =  PBuj  *  PAuj+x  =  Auy+J  for 
j  =  0,...,  k-\.  On  the  other  hand,  (5.14)  for  j  =  k- 1  shows  that  Buk  e  rge  P  -  rge  A  and  thus  that 
Buk  m  Auk+1  for  some  vector  u4+1.  This  completes  the  proof. 

For  the  proof  of  the  equivalence  of  (5.9)  and  (5.10)  we  begin  by  applying  Lemma  5.1  to  (5.9)  which 
corresponds  to  (5.12)  for  A  =  A*_ t,  B  =  Bk. j.  Hence,  we  conclude  that  the  validity  of  (5.9)  is  equivalent  with 
the  condition  that 


A*-i« o  =  0,  Bk.xu0  ~  Ak_xU i,  Bk.tux  m  At_[U i 


(5.15) 
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implies  u0  =  0.  Evidently  (5.15)  is  (5.12)  for  A  =  At_2,  B  =  fl*_2  and  thus  the  lemma  can  be  applied  again.  In 
other  words,  by  repeating  the  process  we  arrive  after  k- 1  applications  of  Lemma  5.1  at  the  condition  (5 AO) 
which  then  completes  the  proof  of  the  equivalence  of  (5.9)  and  (5.10). 

As  before,  let  now  k  be  the  smallest  integer  for  which  (5.9)  holds;  that  is,  for  which  At+1  is  invertible.  If 
no  such  integer  exists  wc  set  k  =  °°.  Then  for  k  <  »  the  system 

&k+\X  —  f  X([q)  —  Xq 

has  a  unique  solution  which,  from  our  theory,  is  a  solution  to  (5.1)  if  and  only  if  A,p0  +  B/xo  = 

j  =  0,  •  •  •  ,  k,  where  p0  is  characterized  by  Ak+lp0  +  Bk+lx0  =  /i+1(fo). 

This  solvability  result  raises,  of  course,  the  question  how  our  condition  (5.9)  (or  (5.10))  relate  to  the  regu¬ 
larity  and  the  index  k  of  the  matrix  pencil  (A ,  B ).  This  is  answered  as  follows: 

Theorem  5.1:  If  the  matrix  pencil  (A ,  B)  is  regular  and  rank  A  <  n  (so  that  K  £  1)  then  k  =  k-1  and  hence 
k  <  oo.  Conversely,  if  k  <  °o  then  (A ,  B )  is  regular  and  hence  k  =  K-1  if  rank  A  <  n . 

Proof:  Let  (A ,  B )  be  regular  and  choose  \  such  that  C  =  B  +  XA  is  invertible.  Then  the  index  k  is  defined  as 

the  smallest  integer  for  which  (5.5)  holds.  Let  uK  be  any  vector  for  which 

(C"lA)’t+1«K  =  0  (5.16) 

and  set  uj  =  C~'Auj +1  =  (C_1A)K^uK,  j  =  k-1,...,0.  Since  (5.16)  implies  that  A(C~lA)KuK  =  0  it  follows  that 
Au0  =  0.  Hence  altogether  we  have  the  equations 

Au0  =  0,  Cuj  =  Au/+1,  j  =  0.....K-1,  (5.17) 

and,  conversely,  (5.16)  holds  whenever  u0,  •  •  •  ,  uK  solves  (5.17).  Now,  by  (5.6),  the  condition  (5.16)  implies 

that  u0  =  (C_lA  )Kut  =  0  for  every  solution  uq,  •  •  •  ,  uK  of  (5.17).  This  implication  is  not  true  when  K  is 

replaced  by  any  smaller  integer  as  is  easily  seen  when  taking  /  <  k  and  choosing  u(  €  ker  C1*1  not  in  iter  C1: 

In  fact,  the  family  Uj  =  C~lAuj+u  j  =  /— 1,  ,0  obviously  satisfies  Au0* 0,  Cuj  •  Au)k.x, 

j  ■  0,  •  •  •  ,  /-l,  but  u0  *  0  since  u  is  not  in  ker  ( C~lA)‘ . 

Thus,  at  this  stage,  we  know  that  k  is  the  smallest  integer  for  which  existence  of  the  solution 
u o,  •  •  •  ,  ux  of  (5.17)  implies  that  u0  =  0.  Given  arbitrary  vectors  uq,  •  •  ■  ,  uK  in  R  ",  define 
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7-> 

V0  =  «0  •  vl  =  Z 

t  =0 


(-X)‘u,_  ,7  =1,  •  •  •  ,  tc 


(5.18) 


With  C  =  B  +  XA  ,  a  straightforward  induction  argument  now  shows  that  if  u0,  •  ■  •  ,  uK  solve  (5.17), 

then 


Av0  =  0  ,  Bvj  =  AvJ+ 1  ,7=0,  •  •  •  ,  K-l  . 


(5.19) 


Conversely,  if  v0,  •  ,  vK  solve  (5.19),  then  the  relation  (5.18)  may  be  inverted  to  produce  a  solution 

uq.  •  •  •  ,  uK  of  (5.17).  In  fact,  we  have 


7-1 

u0  =  V0  ,  Uj  =  £ 
*=0 


7-1 

7 


,7  =  1,  ■  •  •  ,  K 


It  then  clearly  follows  that  k  is  the  smallest  integer  for  which  existence  of  a  solution  v0,  •  •  •  ,  vK  of  (5.19) 
implies  that  v0  =  0,  whence  k  =  k-1  by  definition  of  k. 

Conversely,  suppose  now  that  k  <°°.  Then,  our  theory  ensures  that  (5.1)  has  at  most  one  C2-solution  for 
each  prescribed  initial  condition.  In  particular,  this  holds  for  the  homogeneous  problem 


Ax'  +  Bx  =0,  x(t o)  =  0. 


(5.20) 


But  then  the  pencil  (A,  B)  must  be  regular  because,  otherwise,  the  system  (5.20)  has  infinitely  many  C“ 
solutions  (see  e.g.  [4]). 

The  case  k  =  1  was  noted  already  in  [4],  but  the  general  result  appears  to  be  new. 


6.  Nonlinear  Problems  with  Higher  Index 

As  shown  in  the  previous  section,  for  the  linear  problems  (5.1)  the  sufficient  condition  (2.23)  is  equivalent 
with  the  statement  that  the  pencil  (A ,  B )  is  regular  with  index  one.  The  discussion  in  that  section  suggests  that 
we  may  proceed  analogously  when,  for  the  general  (autonomous)  problem 

F( jt.jO-O,  (6.1) 


the  condition  (2.12)  docs  not  hold  but  F  is  of  class  Cm,m  2  3. 
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The  first  step  in  the  construction  of  a  sequence  of  problems  corresponding  to  (5.9a/b)  was  already  done  in 
Section  2  when  F  was  only  C2.  In  fact,  we  differentiated  (6.1)  and  then  applied  the  projections  P  and  Q  to 
obtain  the  new  problem  (2.9);  that  is, 

G(x,x)  sP(x,  x)F{x,x')  +  Q(x,x')DxF(x,x')x:  =  0.  (6.2) 

Our  sufficient  condition  (2.12)  is  exactly  that  DpG{x,p)  is  invertible  at  the  given  point 
(xo,  Po)  e  F~\ 0)  G~'(0),  and  hence  that,  by  the  implicit  function  theorem,  (6.2)  can  be  transformed  locally 

into  an  explicit  ODE. 

If  this  sufficient  condition  docs  not  hold,  then,  as  in  Section  5,  it  would  now  be  natural  to  set  F°  m  F , 
F 1  ■  G  and,  to  construct  recursively  the  sequence  of  mappings 

FJ+l(x,  p)  =  Pj(x,  p)F'(x,  p)  +  Qj(x,  p)DxFj(x,  p),  j  =  0,  1.  •••  (63) 

where  Pj  again  is  the  orthogonal  projection  onto  rge  DpF‘  and  Qj  =  /„  -  Pj.  Formally,  the  process  is  repeated 
until  the  sufficient  condition  (2.12)  is  satisfied  for  Fk\  that  is,  until  DpFk+l(x ,  p)  is  invertible  at  the  point  under 
consideration. 

As  before,  one  might  consider  calling  the  integer  k+ 1  the  local  index  of  the  problem  at  the  particular 
point.  However,  the  situation  differs  here  in  a  critical  way.  Indeed,  the  very  definition  of  the  iterate  F**1 
assumes  some  smoothness  of  its  predecessor  F‘  and  of  the  projections  Pj  and  Qj.  But  Pj  and  Qj  cannot  even 
be  continuous  at  (x0.  p o)  unless  rank  DpF'(x,p)  is  locally  constant  near  (xo, po)  (in  general,  continuity  of  a 
parametrized  family  in  L(R")  does  not  require  constancy  of  the  rank,  but  it  does  for  projections).  The  validity 
of  such  a  condition  for  j  =  0,  •  •  •  ,  k  is  then  a  necessary  prerequisite  to  iterating  the  procedure  as  outlined 
above.  This  is  not  a  restriction  in  the  linear  case  with  constant  coefficients  of  the  previous  section  because  each 
Fi  involves  matrices  independent  of  (x,  p). 

Conversely,  it  is  easily  seen  that  the  constant  rank  condition  near  (x<>,  p  o)  implies  that  only  one  degree  of 
regularity  is  lost  when  passing  from  F’  to  FJ+\  and  hence  that  F  e  Cm,  m  i  k+ 2,  suffices  to  ensure  that  F*+1 
is  C1.  It  thus  appears  that  aside  from  sufficient  smoothness,  the  constant  rank  condition  near  (xo,  po)  is  the  cru¬ 
cial  ingredient  needed  for  the  definition  of  a  local  index.  Clearly,  in  general,  local  constancy  of  the  rank  cannot 
be  captured  by  a  finite  list  of  requirements  about  F  and  its  derivatives  at  the  point  (x0,  po)  alone:  even  locally. 
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it  is  a  condition  of  a  global  nature. 

Now,  consider  a  problem  with  local  index  k+\  near  (x0,p0).  If  so,  a  C2-solution  to  (6.1)  may  pass 
through  ( x0 ,  p0 )  only  if  F’(x 0,  p0)  =  0,  j  =  0,  •  ■  •  ,  k+ 1.  As  DpFk*l(x0,  po)  e  f«wn  (R*)  by  hypothesis,  the 
problem  Fk+l(x ,  j:')  =  0,  x(r0)  =  *o.  x' o)  -  Po.  may  be  made  into  an  explicit  ODE  for  which  the  solution,  is 
necessarily  C2  if  F  €  Cm ,  m  >  k+  2,  and  is  readily  seen  to  solve  F'(x,  x')  =  0,  j  =  k,  •  •  •  ,  0,  by  recursive 
application  of  Lemma  2.2.  In  sharp  contrast,  nonexistence  of  a  local  index  because  of  failure  of  the  constant 
rank  condition  for  some  iterate  FJ  leads  to  a  singular  ODE  for  which  the  existence  (and  uniqueness)  theory 
should  be  expected  to  be  considerably  different  from  standard  explicit  ODE  theory  in  view  of  the  results  in  [5]. 

In  the  hypothesis  of  the  existence  of  a  global  index;  that  is,  of  the  same  local  index  k+ 1  near  each  point 
of  the  domain  of  definition  E  of  F,  the  global  results  of  Section  3,  corresponding  to  the  case  k  =0,  remain 
valid  since  any  solution  of  (6.1)  is  one  of  Fk(x,x’)  =  0,  a  problem  of  global  index  1.  Finally,  it  should  be 
mentioned  that  problems  of  arbitrary  index  £+1  also  reduce  (locally)  to  explicit  ODE’s  on  manifolds,  at  least 
under  mild  additional  assumptions  such  as  surjectivity  of  the  total  derivatives  DFJ(x,p),  j  -  0,  •  •  •  ,  k.  The 
effect  of  a  higher  index  is  merely  to  shrink  the  dimension  of  the  underlying  manifold:  In  fact,  with  r  -  rank 
DpF(x,p),  r,  =  rank  DpF'(x,p),  j  =  1,  •  •  •  ,k,  we  can  show  that  the  relevant  manifold  has  dimension 
rk  +  ...  +  ri  +  r  -  kn  (compare  with  Section  4  when  k  =  0).  Since  this  dimension  must  be  nonnegative  and 
r,r  i,  •  •  •  ,  rk  <,  n- 1,  one  infers  that,  genetically,  the  index  k+l  cannot  exceed  n.  In  the  linear  case  with  con¬ 
stant  coefficients,  this  result  follows  from  Theorem  5.1  and  k<,  n.  This  example  (via  the  results  in  [4])  also 
shows  that  the  existence  and  uniqueness  theory  for  problems  where  the  iterates  F1  are  defined  beyond  j  =  n 
may  differ  significantly  from  standard  ODE  theory.  A  somewhat  formal  definition  of  the  index  in  terms  of  man¬ 
ifolds  and  related  in  its  spirit  to  the  above  remarks  is  given  in  [6]. 

As  an  example  when  the  constant-rank  condition  remains  valid  for  the  iterated  maps  FJ  we  consider  here 
the  classical  pendulum  problem.  A  pendulum  with  mass  m  attached  at  the  end  of  a  rigid  massless  wire  with 
length  l  attached  at  the  origin  in  the  plane  (xi,  x£  satisfies  the  second-order  DAE 

■*?+*22=/2, 

•  mx\  =  -Xxt  ,  (6A) 

mx2  =  -Xx2  -  mg  , 
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where  X  is  the  (unknown)  tension  of  the  wire  and  g  is  the  gravity  constant. 


In  order  to  reduce  the  problem  to  a  first  order  DAE  we  introduce  the  new  variables  xj  =  and  x4  =  Xj 
and  write  x5  =  k/m.  Then,  the  problem  assumes  the  form  (6.1)  with  x=(xlr  ■■■  ,  xs)eR5, 
p  =  (pi,  •  •  ■  ,  €  R5  and 


F(x,p)  = 


P 1  -  *3 
Pi-** 

Pi  +  x  1X5 
P*  +x2x5  +  g 


(6.5) 


Evidently,  if  ej es  denote  the  standard  basis  vectors  of  Rs  then  rge  DpF (x, p)  =  span  [ej,j  =  2,..,5}. 

Hence  the  conditions  (2.2a/b)  arc  satisfied  and  the  orthogonal  projections  P,Q  onto  rge  DpF(x,p)  and  its 
complement,  respectively,  are  independent  of  ( x,p ).  A  straightforward  calculation  shows  that  the  mapping 
F1  *  G  of  (6.2)  has  the  form 


F\x,p)  = 


2xj>i  +  2X2P2 
P\  ~*3 
IP  2  ~  *4 
P3  +  XjX5 
[P4  +  X2XS  +  g 


(6.6) 


It  is  easily  checked  that  DpFl(x,p )  has  the  constant  rank  4  and  hence  is  not  invertible.  In  other  words,  the 
sufficient  condition  (2.12)  does  not  hold  for  the  pendulum  problem.  But  since  Fl  does  indeed  satisfy  again  the 
constant  rank  condition,  we  may  proceed. 

With  z(x)  =  (1,  -2xi,  -2x2,  0,  0)r  the  orthogonal  projection  onto  [rge  DpFl(x, p)f  is  the  rank-one 
matrix  Qi(x)  =  t{x)z(^)T lz(^)T z{x).  Since  with  Q:  also  Pi  =  1 5  -  Q\  depends  only  on  x  it  follows  that 


DpF2(x, p)u  =Pi(x)DpF'(x,p)«  +  Q\(x)Dp[DxFl(x,p)p]u  .  (6.7) 

Thus  we  have  u  6  ker  DpF2(x,p)  if  and  only  if  both  terms  on  the  right  side  of  (6.7)  are  zero.  Since 
P\DpFx  =DpFx  it  follows  that  dim  ker  DpF2(x,p)  5  dim  leer  DpFl(x,  p)  =  1.  Now  a  short  calculation 
shows  that  DpF\x,  p)e$  =  0  and  therefore  that  rank  DpF\x,  p)  =  4. 

In  order  to  check  whether  the  sufficient  condition  (2.12)  holds  for  F2  that  is,  whether  DpF3fa,p)  maps 
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onto  R5  let  u=ae5  e  ker  DpF2(x ,  p).  Then  D2pFl  ■  0,  and 

DPf[DxFx(x,  p)p](u,  q)  =  (4qlul  +  4q2u2,  0,  0,  0,  0)r  from  the  explicit  calculation  of  DxFx(x ,  p)p  whence, 
because  of  kcr  DpF2{x ,  p)  =  span  {es}  we  have  Dp2pF2(x ,  p)(q ,  e5)  =  0.  Since  £>i(x)e5»0,  we  have 

DQ i(x)«s  =  0  and  DP,{x) c>5  =  -DQl(x)e5  =  0,  Moreover,  one  also  finds  Q\(x)DXfi[DxFl(x,  p)p](p,  es)  =  0. 

Thus  altogether, 

DXfiF2{x,  p)(p,u)  =  P i(x)D2j,F{(x ,  p)(p,u)  =  0  , 

because  DXJ)Fx(x,  p)(p,  =  0  as  a  result  of  DpFl(x,  p)e$  =  0.  Similarly,  we  have 

DxF\x,  p)u  =  Pl(x)DxF\x,  p)u  =  DxF\x,  p)u  because  DxFl(x ,  p)e$  =  ( 0,  0,  0,  x1(  x?)T  ■  v(x)  and  hence 
Qi(x)v(x)  =  0. 

Thus  in  this  case  the  sufficient  condition  (2.12)  for  F2  simply  requires  that  when 

ctv(x)  €  rge  DpF2(x,p)  then  a  =  0.  Since  the  two  terms  on  the  right  of  (6.7)  are  complementary  and 
Qt(x)v(x)  =  0  this  means  that  the  equations  DpFl(x,  p)q  =  av(x)  and  Q\{x)Dp[DxFl(x ,  p)p]q  =  0  only  have 
a  solution  q  when  a  =  0.  In  fact,  the  first  of  these  equations  amounts  to  q  j  =  q2  =  0,  q3  =  CLr1(  and  q4  =  ax2, 
while  from  the  second  equations  it  follows  that  q 3  =  qA  =  0.  Hence,  for  (xlt  xj  *  (0,  0)  we  indeed  have  a  =  0 
as  desired. 

Because  F 1  and  F2  satisfy  the  constant-rank  condition,  the  existence  and  uniqueness  results  apply  for 
every  initial  condition  (x0,  p0)e  R5  satisfying  FJ(x0,  p0)  =  0,  j  =  0,  •  •  •  ,3  since  (x01,  X02)  *  (0,  0)  from 
F(xo<  Po)  =  0  (see  (6.5)).  Moreover,  here  we  may  indeed  say  that  the  problem  has  global  index  three  since  the 
condition  (x  1,  xj)  *  (0,  0)  is  not  a  restriction  along  and  hence  near  the  solutions  of  (6.4). 
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